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We report on an interference effect in multiple scattering by resonant scatterers resulting in 
enhanced forward scattering, violating Ohm's law for photons. The underlying mechanism of this 
wave effect is superradiance, which we have investigated using cold atoms as a toy model. We 
present numerical and experimental evidences for this superradiant forward scattering, which is 
robust against disorder and configuration averaging. 

PACS numbers: 42.50.Gy, 42.25.Bs, 03.65.Yz , 32.80.Qk 



A rigorous investigation of light propagation in the 
presence of many scatterers requires taking into account 
the effects of interferences. Fortunately in most situa- 
tions interference effects can be neglected and a radia- 
tive transfer equation is thus often used in optics, allow- 
ing a practical approach to scattering of light in com- 
plex media. Some exceptions are however known, where 
interferences are fundamental for the observed features. 
One prominent feature is coherent backscattering (CBS) 
^ [5] . Interferences also come into play in experiments 
on universal conductance fluctuations |3j and more recent 
experiments probing the local density of states [IHS] . A 
very active line of research concerns strong localization 
of waves of light [7J [5] ■ However, beyond these examples 
few features where interference effects survive after con- 
figuration average in optical dense media are known. In 
the limit of very dense samples, close to the strong local- 
ization regime, another hitherto non observed features is 
coherent forward scattering [9j. Most of the above fea- 
tures are genuine wave effects and can be present in the 
propagation of different types of waves, such as matter 
waves [TUIIT3] or acoustics [T4"1[T5]. 

If scattering of light is enhanced by a narrow resonance, 
additional features appear. Indeed at resonance the po- 
larisability becomes imaginary, describing the opening of 
the system to many vacuum modes. As a result cooper- 
ative scattering as pioneered by Dicke [H] leads to su- 
perradiance [T7HT9] l23l [24] investigated in various recent 
experiments J2QH22J . Analytical investigations are diffi- 
cult, as a rigorous theory for wave propagation in mul- 
tiple scattering is not developed. In order to go beyond 
numerical studies, one either resorts to an effective index 
approach [25] or a diagrammatic analysis [26 . For reso- 
nant scattering, these approaches are however not valid 
[2"7] and the situation is reminiscent of resonant scatter- 
ing by matter waves at the unitary limit [35] [22] ■ 

In this letter we present a new interference effect in 
multiple scattering. We start our investigation of this 
steady state superradiant forward scattering (SFS) in 
the multiple scattering regime using numerical simula- 
tions before presenting experimental results of such phe- 



nomenon. 

The relevance of interference effects in multiple scatter- 
ing is best highlighted by a comparison to an incoherent 
model, where the propagation of photons are described 
by classical ray optics. We investigate the role of inter- 
ferences using different wave models. The first one is a 
microscopic model based on N coupled dipole equations 
[30) . a type of many body solution for our wave equa- 
tion. We also use two analytical models which amount 
to a rough mean field Ansatz: the driven timed Dicke 
model [TBI I3T] , where the N randomly-distributed parti- 
cles have the same amplitude and a phase imprinted by 
the laser, successfully used to describe previous experi- 
ments and a partial wave approach [35] used to predict 
Mic resonances. We choose to use a simplified scalar 
model (instead of the more rigorous vectorial model [T7| ) 
as this highlights the essence of SFS which does not rely 
on the specific vectorial nature of light. 

A textbook situation to study the propagation of light 
through disordered media is the slab configuration. We 
focus our interest in this letter on resonant scattering 
and define the optical thickness bo — ncrgL, where n 
is the spatial density of atoms, <7o = 3A 2 /27r the res- 
onant cross section at the wavelength A for individual 
two-level atoms and kL the thickness of the slab in units 
of /c _1 = A/2-7T. In the absence of cooperative effects, 
this optical thickness is related to the forward (ballistic) 
transmission via T co h = e~ b ° . We compute the intensity 
distribution inside the slab and the diffuse transmission 
through the sample using different numerical methods. 
As one can see in Fig. [TJ the intensity distribution ob- 
tained from our many body solution (blue squares in Fig. 
[TJ inside the sample for an optical thickness of bo = 8 is 
far from being uniform. This might a priori not come as 
a surprise, as an incoherent model of multiple scattering 
of photons undergoing a random walk inside the sample 
provides a similar intensity distribution (black circles in 
Fig. [TJ. For this incoherent model we perform a Monte 
Carlo simulation, isotropically scattering photons after 
a scattering mean free path Z sc = L/bo, where L is the 
thickness of the slab. We note that for the many body 
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FIG. 1: (color online). Intensity distribution inside a slab 
with thickness of atoms with a thickness L z = 20/fc and bo = 8 
and a transverse size of L x = L y = 80/fc. The local inten- 
sity inside the slab is computed using an incoherent model 
(black spheres) and compared to the excited state population 
derived in a driven timed Dicke model (green triangles) and a 
partial wave Ansatz (pink diamonds). The limit of indepen- 
dent dipoles (yellow stars) is recovered in the driven timed 
Dicke model for f>o — t 0. The many body solution is obtained 
using N coupled dipole equations (blue squares). 

solution, the local field has important fluctuations reflect- 
ing the speckle nature of the field in multiple scattering. 
For the data shown in Fig. [T] we have thus used an in- 
tegration over the transverse directions inside the slab, 
applied a coarse graining with a resolution of k5z — 2 
and averaged over 30 configurations. As one can see in 
Fig. [T] our driven timed Dicke model very poorly de- 
scribes the intensity distribution inside the sample: for 
the slab geometry considered here, the population of the 
excited state (green triangles in Fig. fl]), homogeneous 
in this Ansatz, scales as (1 + &o / 6) — 2 33J . When using a 
dilute sample (vanishing 6 ) this model recovers the in- 
dependent atom case (yellow stars in Fig. [I]). For large 
&o however only the first layer of the sample is excited 
by the incident field alone. Dipoles located in the deeper 
layers of the sample are excited by the interference of the 
incident field with the field scattered by all other dipoles 
and the homogeneous Ansatz for dipole excitation is no 
longer appropriate. We also show the result of a partial 
wave approach (pink diamonds in Fig. [IJ, expected to be 
valid not only when the dipole are excited by the exter- 
nal field. However this model assumes an effective index 
response. The imaginary part of the index of refraction 
does not discriminate between absorption and scatter- 
ing and only contributes to an extinction of the incident 
wave, according to the Beer-Lambert law. Whereas we 
expect the use of an effective index of refraction to be well 
suited for off-resonant scattering, where the atomic po- 
larisability is almost real, such an effective field model has 
severe limitations in the presence of resonant scattering. 



Similar situations are also encountered in metamaterials, 
with different mechanisms at the origin of the resonances 

Let us now analyze a better known quantity in trans- 
port through random media: the diffuse transmission 
Tdiff through a disordered sample. In Fig. ([2]), we show 
the result for T^is using the two models which seemed 
to provide similar results for the intensity distribution 
inside the sample: the incoherent model and our many 
body simulation. For large optical thickness, the diffuse 
transmission is described in the incoherent model by a 
diffusion equation reminiscent of the Drude model for 
electrons transmitted through a wire. As a result this 
model predicts for large bo'. T^m oc l/6o, dubbed Ohm's 
law for photons (black spheres in Fig. [2j. We note that 
energy conservation insures that T^g + i?diff + T co b = 1, 
where i?diff is the diffuse reflection coefficient and T co h the 
remaining unscattered coherent forward beam. In con- 
trast to the result of the incoherent model, amounting to 
neglect all interference effects in multiple scattering, our 
many body solution provides striking different results for 
Tdiff (blue squares in Fig. [2]). We have normalized the 
scattered light across the forward boundary of the slab to 
the total scattered power (which approaches a constant 
value at large bO) and checked for energy conservation 
at low optical thickness. The surprising result is that, 
instead of a decrease, we observe a saturation for T di ff 
for large bo, as if the "resistance" for photons to cross 
this sample would vanish. Resonant scattering of light 
through a large sample thus violates Ohm's law for pho- 
tons. We note that in addition to this saturation of the 
diffuse transmission at a non zero value, the many body 
simulation also predicts a different slope at small bo- We 
attribute this feature, which can be assimilated to a re- 
duced effective optical thickness in the many body simu- 
lation, to the consequence of a long range dipole blockade 

Further insight can be obtained from the emission di- 
agram. As the cold atom sample used in our experimen- 
tal setup has a spherical symmetry with a radial Gaus- 
sian density distribution n(r) oc e~ r /( 2<t r) and a opti- 
cal thickness of bo = J n(0, 0, z)aodz = 3N/(ka r ) 2 we 
compute the emission diagram for this specific geome- 
try in Fig. [3j The incoherent model (black line in Fig. 
[3]), where photons are isotropically scattered after a dis- 
tance l/(n(r)<7o), yields a smooth intensity profile, with 
slightly more photons scattered backwards than forward. 
We note that for this Gaussian sphere, many photons 
can be scattered at large distance from the center where 
the optical thickness is small, yielding a less pronounced 
backward/forward ratio than a slab geometry. The result 
of the driven timed Dicke model (green line in Fig. [3]) in- 
cludes the pronounced forward emission lobe as expected 
from the cooperative interference effect of the scattered 
light. The more rigorous result of the many body solu- 
tion (blue line in Fig. |3| in addition includes CBS. This 
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FIG. 2: (color online). Ohm's law for photons. Diffuse trans- 
mission of a plane wave through a slab of increasing optical 
thickness bo, for the incoherent model (black spheres) and the 
many body solution (blue squares). 



feature, observed in experiments more than 10 years ago 
[23 [36] is not described by the driven timed Dicke model. 
Indeed CBS requires dipoles to be excited by the field 
scattered by others dipoles and not by the incident field 
alone. In contrast, diagrammatic models which describe 
CBS by atomic scatterers [ST] do not include the physics 
of superradiance. SFS is thus beyond the reach of such 
models used in the past to describe resonant multiple 
scattering. 

Possible experimental signatures of SFS obviously in- 
clude the large diffuse transmission around forward direc- 
tion. Even though enhanced scattering around backward 
direction (CBS) has been observed in experiments by di- 
rect detection of the scattered light, the direct detection 
of SFS around forward direction is more challening. In- 
deed, despite that fact that SFS can be many orders of 
magnitude larger than CBS, stray scattering by imper- 
fections in experiments often makes direct light detection 
around forward direction more difficult. We therefore 
chose to explore a different signature of SFS, exploit- 
ing the potential of cold atoms to be an efficient detec- 
tor of light scattering [2D]. As one can anticipate from 
Fig. [3] incoherent scattering would increase the average 
radiation pressure force, as more photons are scattered 
in backward direction. SFS on the contrary leads to a 
drastic reduction of the radiation pressure force at reso- 
nance, as many photons are scattered around the forward 
direction. We note that CBS also contributes to a mod- 
ification of the average radiation pressure force (slightly 
increasing the force), but the amount of light scattered in 
the CBS cone is so small that its detection via a change 
of the radiation pressure would be a tour de force. 

We thus apply the following experimental procedure 
to probe SFS in cold atoms. First, we load a magneto- 
optical trap (MOT) with 3 x 10 7 atoms of 87 Rb in 50 
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FIG. 3: (color online). Emission diagram obtained numer- 
ically for a Gaussian sphere (k<jR = 20) of central optical 
thickness bo = 10 for the incoherent model using 13000 pho- 
tons (black line), the timed Dicke Ansatz (green line) and the 
many body simulation (blue line), averaged over 200 configu- 
rations. Logarithmic scale is used for the intensity. 



ms using the setup described in [20 . We then apply a 
50 ms temporal dark MOT period where the intensity of 
the repumping laser is reduced by a factor of 10 and the 
cooling laser is tuned to — 10 T from the F = 2 —> F' = 3 
D2 line. This allows to compress the cloud and to pro- 
duce a smooth Gaussian shaped distribution of atoms. 
To control the optical thickness at the end of this dark 
MOT period, the repumper detuning is varied between 
— 7F to — 2.5 r, keeping the desired amount of atoms in 
the F — 2 state without affecting the size (ctr — 270 /rat, 
kaji w 2 x 10 3 ), shape and temperature (~ 20 fiK) of 
the cloud. We then switch off all laser beams and mag- 
netic field gradients, leaving the atoms in free fall. We 
then apply an horizontal, circularly-polarized "pushing" 
beam, tuned close to the F = 2 — >• F' = 3 transition 
for 50 fis. The pushing beam has a waist wq — 12 mm 
and its carefully calibrated intensity is adjusted to have a 
saturation parameter s = 10~ 2 . Each atom in the F = 2 
state scatters on average 10 photons. Such a small num- 
ber of scattered photons prevents from any depumping 
effect into the F = 1 state during the pushing process. 
After a time of flight expansion of 12 ms, we image the 
position of the atomic cloud via standard off resonant 
(detuning A s» — 2T) absorption scheme. A Gaussian fit 
of the absorption image yields the position of the center 
of mass of the atomic cloud after time of flight and thus 
the average radiation pressure force. Each experimental 
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FIG. 4: (color online). Force ratio as the function of optical 
density bo (log-log scale) for resonant laser light. 



point (see figure[4]) is an average over 10 realizations with 
error-bars reflecting the shot to shot fluctuations of the 
experiment. We normalize the measured average radia- 
tion pressure force by the single atom force -Fi , where F\ 
is computed without any adjustable parameter from the 
measured intensity of the pushing which (known within 
5 % accuracy and taking into account losses by the vac- 
uum windows). The experimental value of the optical 
thickness bo is obtained by standard absorption imaging. 

In Fig. [4] we show the signature of SFS on the aver- 
age radiation pressure force. As shown by the numerical 
results of our many body solution (blue squares), the 
net result of SFS on the average radiation pressure is a 
strong reduction of the force. Indeed, as most photons 
are crossing the sample close to forward direction, the 
photons leave the system without important transfer of 
momentum. The larger bo, the more photons are scat- 
tering in the forward lobe of width A8 « l/(kaji). The 
experimental signature of our study of SFS (red spheres 
with statistical errors bars) clearly shows an important 
reduction of the average radiation pressure force. This 
is in striking contrast to the prediction of the incoherent 
model (black spheres). Indeed, if one could neglect all in- 
terference effects for the propagation of light through this 
optically thick sample of atoms, more and more photons 
would be scattered backward. For very large bo, almost 
no photon would cross the sample (bare a few on the 
edges of the sample) and the average radiation pressure 
force would be above that of the independent atom limit. 
We also indicate on Fig. [4]the prediction using our driven 
timed Dicke model. We stress that despite the apparent 
agreement with the experimental values, we do not con- 
sider this model to provide an accurate account of the 
detailed mechanism at work for SFS. Indeed, as shown 
in Fig. [T] this mean field result does not describe cor- 
rectly the intensity distribution inside the sample. 

We note that despite the agreement between the exper- 



imental results and the many body solution, one needs to 
keep in mind that small systematic corrections could be 
expected for the calibration of our values of bo- Indeed 
we have used a standard low intensity absorption tech- 
nique to measure bo . This does not take into account 
the differences between the transition strengths for var- 
ious Zeeman sublevels, which could be taken care of by 
more high intensity absorption imaging [38 . Also, as 
discussed above, our many body solution is an approx- 
imate scalar model. Even if we expect this to capture 
the main features, in particular for dilute cloud of atoms 
where near field effects can be neglected, the emission di- 
agram will have some polarization dependence not taken 
into account by any of the models presented in this work. 
Minor additional effects which could affect SFS is reduced 
spatial or temporal coherence as e.g. induced by finite 
additional disorder, temperature or saturation [24] , 

In summary, we have introduced an interference ef- 
fect in multiple scattering based on superradiance. In 
SFS, Ohm's law for photons is violated and most pho- 
tons are transmitted through the sample, as if cooper- 
ativity would turn a diffuse medium (with a resistance 
for the transmission of photons) into a superconductor 
of light. We speculate that similar resonant effects could 
be induced with different type of resonances (internal res- 
onances, shape resonances or confinement-induced reso- 
nances) and applied to different waves (ultra-cold atoms, 
electrons) . We also note that the long range dipole-dipole 
coupling used in our many body simulations bears some 
resemblance with the Heisenberg model used at the ba- 
sis of the resonating valence bond theory [35] and that 
light scattering by cold atoms might provide a useful toy- 
model to further study these effects. 
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